THE 1/2 COMPLEX BRUNO FUNCTION AND THE YOCCOZ 
FUNCTION. A NUMERICAL STUDY OF THE 
MARMI MOUSSA YOCCOZ CONJECTURE. 



TIMOTEO CARLETTI 

Abstract. We study the 1/2— Complex Bruno function and we produce an 
algorithm to evaluate it numerically, giving a characterization of the monoid 
M = Mt U Ms- We use this algorithm to test the Marmi-Moussa-Yoccoz 
Conjecture about the Holder continuity of the function z \— > — i B(-z)+log C/(e 27rIZ ) 
on {z £ C : ^5z > 0}, where B is the 1/2— complex Bruno function and U is 
the Yoccoz function. We give a positive answer to an explicit question of S. 
Marmi et al IMarmi et al. 20011 . 



1. Introduction 

The real Bruno functions are arithmetical functions B a : R\Q — > R+ U{+oo}, 
a 6 [1/2, 1] which characterize numbers by their rate of approximation by rationals. 
They have been introduced by J.-C. Yoccoz [Yoccoz 1995j (cases a = 1/2 and 
a = 1) and then studied in a more general context in jMarmi et al. 1997) . 

For their relationship with arithmetical properties of real numbers, Bruno's func- 
tions enter in a huge number of dynamical system problems involving small divi- 
sors, for instance in the problem of the stability of a fixed point of a holomor- 
phic diffeomorphism of a complex variable (the so called Schrdder-Siegel prob- 
lem) | Yoccoz 1995] . in the Schrdder-Siegel problem in the Gevrey setting in one 
complex variable [Carletti-Marmi 2f)f)0| and in some local conjugacy problems: 
Semistandard map M armi 1990llDavie 1994| . analytic circle diffeomorphisms |Yoccoz 2002 
and some analytic area-preserving annulus map including the Standard map and 
some of its generalizations Berretti Gen tile 2001| . 

Let us now concentrate on the 1/2-Bruno function . -B1/2 is Z-periodic, even 
(for this reason it is also called even Bruno 's function) and verifies the functional 
equation: 

(1.1) B 1/2 (x) = -logx + xB 1/2 (x- 1 ) ze (0,1/2). 

The set B = {x € R : B 1 n (x) < +00} is called the set of Bruno's numbers: by (|1.1|) 
it follows that B is invariant under the action of the modular group GL(2,Z) = 

{ (c d) ' a )b,c,d € Z, ad — be = ±lj- The Bruno function can be extended to 

rational numbers by setting Bi/ 2 (x) — +00 when x G Q. 

Date: February 1, 2008. 

Key words and phrases. Complex Bruno Function, Yoccoz Function, linearization of quadratic 
polynomial, Littlewood Paley dyadic decomposition, continued fraction, Farey series. 

^From IMarmi et al. 19971 we know that the difference of any two Bruno's functions is in 
i°°(M). 



2 



TIMOTEO CARLETTI 



Using the continued fraction algorithm one can solve to obtain: 
(1-2) B 1/2 (x) = (x) log x' k \ 

k>0 

where xo — x, Xk = A l / 2 (%k-i), P-i = 1, ftk = W.j—Q x i an d A X j 2 is the nearest 
integer Continued Fraction map. In S l2.1l we will give a brief account of useful facts 
concerning continued fractions. 

In [Marmi et al. 2001| the complex Bruno function has been introduced 2 , more 
precisely authors defined an analytic map B : HI + — > H + , where H + is the upper 
Poincare half plane, Z~periodic, which verifies a functional equation similar to the 
one for the 1-Bruno function. The boundary behavior of B is given by (see Theorem 
5.19 and §5.2.9 of |Marmi et al. 2001| 1: 

(1) let H > 0, then the imaginary part of H(z) tends to Bi/ 2 (3iz) when — > 
and z G {C G H+ : 3C > \^C\ H }, whenever 5iz G B; 

(2) 9?B(z) is bounded on H + , its trace is continuous at irrational points and 
it has a jump of ir/q for = p/q G Q. 

In section [3] we introduce an explicit formula for the 1 /2-complex Bruno function 
which corrects a small error in §A.4.4 page 836, and gives more details than Ap- 
pendix A. 4 of Ma rmi et al. 200T] . We will also give an algorithm to compute it 
numerically. 

1.1. The Yoccoz function. We already observed that the function Bi/ 2 is related 
to the stability problem of a fixed point of an analytic diffeomorphism of C, in 
the rest of this section we will show this relation by describing the Yoccoz result 
( [Yoccoz 1995| . Chapter II). Let A G C* and let us consider the quadratic polynomial 
P\(z) = Xz(l — z). The origin is a fixed point and we are interested in studying its 
stability. If |A| < 1 (hyperbolic case) then it follows from the results of Poincare 
and Koenigs that the origin is stable, whereas if A = e 2j!lp / q (parabolic case) the 
origin is not stable. 

Let now consider A G D* and let H x (z) be the conformal map which locally 
linearizes P\ (its existence is guaranteed by the Poincare-Koenigs results) : 

(1.3) P x oH x = H x oR x , 

where R x {z) — Xz, and let us denote by ^(A) the radius of convergence of H x . 

One can prove that H x can be analytically continued to a larger set, the basin 
of attraction of 0: {z G C : P x n {z) — > 0,n — > +oo}, but not to the whole of C, 
and it has a unique singular point on its circle of convergence D r2 (A)7 which will be 
denoted by U(X) G C- The function U : D* — > C is called the Yoccoz function. 

Yoccoz proved that U has an analytic bounded extension to D and moreover it 
can be obtained as limit of polynomials U n (X) = X~ n P£ n (z cr it), uniformly over 
compact subsets of B, where z cr it = 1/2 is the critical point of the quadratic 
polynomial. Being this extension not identically zero, by a classical result of Fatou, 
the Yoccoz function has radial limits almost everywhere, and the set Ao G S 1 for 
which limsup A ^ Ao U(X) — has zero measure. This means that the quadratic 

2 Following the notation introduced for the real Bruno functions, we should call this complex 
extension the 1-complex Bruno function. In fact we will see at the end of section |5] that it 
is constructed "following" the Gauss continued fraction algorithm. In this way we could also 
distinguish it from the 1 /2-complex Bruno function that we will introduce in sectional "following" 
the nearest integer continued fraction algorithm. 
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polynomial is linearizable (r2(Ao) = |C/(Ao)| > 0) for a full measure set of Ao £ S 1 , 
but the proof doesn't give any information on this set. 

When |A| = 1 and A is not a root of the unity, assuming A = e 27rluJ , for some 
irrational \u>\ < 1/2, Yoccoz proved f [Yoccoz 1995| . Theorem 1.8 Chapter II) that 
P\(z) is linearizable if and only if uo £ moreover there exists a constant C\, 

and for all e > a constant C(e) such that for all u> £ Bi/ 2 : 

Ci < logr 2 (e 2 ^) + B 1/2 (uo) < C(e) + eB 1/2 (co) . 

We are then interested in studying the function uo i— > log \ U (e 27 ™") | + B 1 / 2 (uo) and 
some "natural" questions arise ( \ Yoccoz 1995) §3.2 page 72): 

(Yoccoz's Conjecture). Is the function id i— > log \U (e 2muj ) \+B 1 / 2 (uo) 

bounded for 

Motivated by numerical results of |Marmi 1990| and by some analytic properties of 
the real Bruno function (see the following Remark 1 1.1 1 and [Marmi et al. 1997 | it 
has been conjectured that: 

(Marmi-Moussa- Yoccoz's Conjecture). The function, defined on 
the set of Bruno number, uo i— ► log | C/ (e 27 ™") | + Bi/ 2 (uo), extends 
to a 1/2-Holder continuous function on R. 

Remark 1.1 (Why l/2-H61der?). In Mar mi et al. 1997] authors proved a "sta- 
bility result" for B\j 2 (§4 page 285). Let us rewrite the functional equation for the 
\ 1 2-Bruno function as follow: 

[B 1/2 (x) - xB 1/2 (se -1 )] =-loga;, 

if we add to the r.h.s. a "regular term" f , say rj-Holder continuous, and we call 
Bf the solution of: 

[B f (x) - xB f (x- 1 )] = - logx + f{x) , 
then Bij 2 — Bf is 1/2-Holder continuous if f is at least 1/2-Holder. Hence if 



we prove 3 that the function u> 



log uj 



\og\U(e 2mu ) |- a; log U (e 2 ™ 

is Holder continuous with exponent r\ > 1/2, for uo £ [0, 1/2], then the Conjecture 
holds. 



Very recently X. Buff and A. Cheritat jBjifcCheritat 2003) proved the Yoccoz 
conjecture. Whereas the Marmi-Moussa- Yoccoz conjecture is still open. We will 
be interested in the following conjecture, equivalent to the one of Marmi-Moussa- 
Yoccoz: 

(Marmi-Moussa- Yoccoz's Conjecture) The analytic function, de- 
fined on the upper Poincare half plane, z i— ► TL{z) — log U(e 27rtz ) — 
i B(z), extends to a 1/2-Holder continuous function on H + . 

The aim of this paper is twofold. First give more insight into the 1/2-complex 
Bruno function and second to make a first step toward the understanding of the last 
conjecture. Our numerical results allows to conclude that TL is ry-Holder continuous 
and we obtain an estimate of the Holder exponent 77 = 0.498 ± 0.004. This gives 
us good numerical evidence that the Marmi-Moussa- Yoccoz conjecture should be 
true. 



■^Transform a function according to ip(x) t—> ip(x) — xip(l/x) to "reduce the strength of singu- 
larities" is the main idea of the Modular Smoothing. We refer to Bu ric et al. 1990l where authors 
describe the method and apply it to the critical function of the Semistandard Map. 
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The paper is organized as follows: in section [3 we introduce the 1/2-complex 
Bruno function and some results from number theory (approximations of rationals 
by rationals), to obtain an algorithm to compute the complex Bruno function. In 
section |21 we explain how to calculate the Yoccoz function and then, after a brief 
introduction of the Littlewood-Paley Theory in § 0J used to test the Holder conti- 
nuity, we present our results in section [J] Appendix 1X1 collects some considerations 
related to technical aspects of our numerical test. 

Acknowledgments. I am grateful to Jacques Laskar for putting at my disposal a 
large amount of CPU times on computers of Astronomie et Systemes Dynamiques 
Team at IMCCE Paris. 



2. The 1/2-Complex Bruno functions 

The aim of this section is to introduce, starting from Appendix A. 4 of [Marmi et al. 2001| 
a complex extension of the 1/2-real Bruno function and to give an algorithm to 
compute it numerically. 

Let us consider / £ L 2 ([0, 1/2]), extended: 1-periodic, f(x + 1) = f(x) for all 
and even f(x) = f(~x) for all x € [—1/2,0], and then let us introduce the 
operator T acting on such / by: 

(2.1) Tf(x)=xf(- 



x / 

we remark that the functional equation (|1.1|) can be rewritten as: 
(2.2) (l-T)B 1/2 (x) = -logic VxG (0,1/2). 

Let (T m ) m >2 be the operators defined by: 



(2.3) (T m f) (x) 




m' m-1/2 

otherwise, 



branch m" 1 
branch to" 



then using the periodicity and the evenness of /, we can rewrite l|2.1|) as follows: 
(2.4) Tf(x) = {*/ (~ - m ) + *f ( m + 1 - ~) } • 



m>2 



To introduce the 1/2-complex Bruno function we have to extend l|2. 4f> to com- 
plex analytic functions; this is done |Marmi et al. 2001] by considering the com- 
plex vector space of holomorphic functions in C \ [0,1/2], vanishing at infinity: 
O 1 (C \ [0, 1/2]) (which is isomorphic to the space of hyperfunctions with support 
contained in [0, 1/2]). So, let ip be the Hilbcrt transform of /: 

( \ 1 f 1/2 f& A 

Viz) = ~ / 



7T Jo X — Z 

then starting from H2.4J) we define the action of T on ip as follows: 
(2.5) T<p(z)=J2 L 9(m)(l + L<r)v(z), 



m>2 
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where g(m) = ( ° £ ), a = ( - 1 \ ) and L^ a ^ acts on O 1 (C \ [0, 1/2]) by: 



(2-6) b^(z) = (a- cz) 



( dz — b 



zd — 6c 



\ a — cz 

In the spirit of (|2.2II we want to consider (1 — T) 1 acting on some ip G O 1 ((£ \ [0, 1/2]), 
and to obtain a Z-periodic, "even function" 4 , we will consider: 

s-l" 



(2.7) 



1 + _L CT ) (1 — T) cp(z-n). 



Let us introduce the operator f defined by (1 + L a )T = T (1 + L CT ) , then from (|2.5|) 

and the relation: (1 — T) 1 = X)»->o T r , we can expand ^1 — f^j in terms of ma- 
trices <?(m) and a, to obtain a sum of matrices of the form: eog(mi) . . . e r ^ig(m r ), 
where r > 1, rrii > 2 and e,_i G {1, cr}, for 1 < i < r. 
Let us set = {1} and for r > 1: 

7W (r) = {.9 G G£(2, Z) : 3e , . . . , e r _i G {1, cr}, mi, . . . , m r > 2 : 

(2.8) g = e j(mi) .. . e r _ig(m r ) |, 

and finally = U r > .M (r) : the 1/2-Monoid (we left to §0 a more detailed 
discussion of this monoid and the reason of its name). 

It remains to specify the "good" <p G O 1 (C \ [0, 1/2]) to apply (|2.7|l . to have the 
wanted properties for B- This is done by considering the Hilbert transform of the 
even, 1-periodic logarithm function defined in (0, 1/2], namely: 



Vl/2(*) 



(2.9) 



1 


,1/2 


— logo; 


rf.r 






7T 


J() 


£ — Z 










1 




— Li 2 




1— 




7T 




7T 



1 

2z 



where Liz (z) is the dilogarithm function |Oesterle 1993) : the analytic continuation 



function to be: 
(2.10) 



2 , to C\[l, +oo). We are now able to define the 1/2-complex Bruno 



B(*) = E 



931/2(2: - n) . 



E Mi+£*) 

geM 

This formula defines 5 an holomorphic function, defined in H + , Z-periodic, with 
an even imaginary part for -Rz G [0, 1/2], and then 6 an odd real part, on the same 
domain. 

Remark 2.1. The 1/2-complex Bruno function is 1-periodic and so we can con- 
sider its Fourier series: B (z) = X^ez he 2mlz ■ Introducing the variable w = e 
the Bruno function is mapped into an analytic function, B (w), defined in D 
which can be extended by continuity to p. Its Taylor series at the origin is B (w) 



27tiz 



4 Here and in the following by even complex function we will mean even w.r.t. SRz — » — Sftz. 

''This claim can be obtained by slightly modification of the proof given in IMarmi et al. 20011 
for the 1-Complex Bruno function and we omit it referring to [Marini ct al. 2001 for any details. 

^This is a standard result for harmonic conjugate functions. For this reason we prefer to speak 
of 1/2-Complex Bruno's function instead of even Complex Bruno's function. 
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X^zeN ' w '> hence Fourier coefficients ofJi(z) corresponding to negative modes are 
all identically zero. Moreover, because of the parity properties o/SftB and SB, its 
Fourier coefficients are all purely imaginary, in fact: 
/■V 2 

k = 2i / [- sin (2irlx) K B (x + it) + cos (2nlx) 3B(i + it)] dx . 
Jo 

The goal of the next sections will be to express <|2.1U[) in terms of a sum over 
a class of rational numbers in such a way we could give (§ 12.41) an algorithm to 
compute it. This will be accomplished thanks to a new characterization (§ 12. 3|) of 
the 1/2-Monoid M., after having introduced some results from Number Theory: 
§ 12.11 Continued Fraction Theory and S I2.2I Farev Series. 

2.1. Continued Fraction. We consider the so called nearest integer Continued 
Fraction algorithm 7 . We state here some basic facts we will need in the follow- 
ing and we refer to |Hardy-Wright 1979| IMarmi et al. 1997] for a more complete 
discussion. Let ||a;|| = min p6Z {x < 1/2 +p}, then to each x £ R we associate a 
continued fraction as follows: 



(2.11) oo = ||x||, 3c = |a;-ao|, e = 



+1 iff x > ao 
— 1 otherwise 



and then inductively for all n > 0, as long as x n ^ 0: 
(2.12) 

-1 iff x^ 1 > a n+ i 



Cri+l — \ x n ~ a n+l\ = A\j2{Xn) , E n+ i — 



-1 otherwise 



We will use the standard compact notation to denote the continued fraction x — 
[(ao,£o) ; ■ ■ ■ , (a« + £nX n ,£n)]- From the definition it follows that x n > 2 and so 
a n > 2. 

Remark 2.2 (Standard form for finite continued fraction). Let [(ao, eo), . . . , (an, En)] 
be a finite continued fraction of length n. Then, whenever a fi = 2, we must also 
have En— i = +1, namely [(ao, Eq), . . . , (a n -i, — 1), (2, +1)] represents the same ra- 
tional number that [(ao, Eo), ■ ■ ■ , (a n —\ — 1, +1), (2, +1)] . Moreover a finite continued 
fraction cannot contain a couple (ai,Ei) — (2, —1) for any I < n. 

We recall, without proof some known results: 

• the continued fraction algorithm stops if and only if x € K \ Q, this corre- 
spondence in bijective up to the standard convention of Remark 12. 21 

• For any positive integer n (or smaller than the length of the finite continued 
fraction) the n th convergent is defined by: 

(2.13) — = [(ao,£ ), ■ ■ ■ , (a n ,E n )\ , 

q?i 



(2.14) 



one can prove that p n and q n are recursively defined by: 

\Pn = CLnPn-1 +£n-lp7i-2 

[<7n = a n q n -i + e n -xq n -2 , 
starting with p_i = o_ 2 = 1, p_ 2 = q_\ = and e_j = 1. 



7 In IMarmi et al. 19971 a one parameter family of continued fraction developments has been 
introduced. The nearest integer continued fraction corresponds to the value 1/2 of the parameter, 
so we will also call it 1 /2-continued fraction. 
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• for all n we have: q n p n -i - p n q n -i = (-l) n £o • ■ ■ £n-i- 

2.2. The Farey Series. Let n € N*, the Farey Series |Hardy-Wright 1979| of 
order n is the set of irreducible fractions in [0, 1] whose denominators do not exceed 
n 8 : 

(2.15) T n = {p/q £ [0, 1] : (p, q) = 1 and q < n}. 

The cardinality of T n is given by the Euler 4>{n) function and so it is asymptotic 
to 3n 2 /7T 2 for n large. The Farey Series is characterized by the following two 
properties |Hardy- Wright 1979| : 

Theorem 2.3. Let n > 1. If p/q and p'/q' are two successive elements of T n , 
then: qp' — q'p = 1 . 

Theorem 2.4. Let n > 1. If p'/q', p/q and p"/q" are three successive elements 
(in this order) of T n , then: 

p _ p' + p" 
q~ q' + q" ' 

Using an idea contained in the proof of Theorem l2 . 41 given in IHardy-Wright 1979| , 
we construct an algorithm (easily implementable on a computer) which allows us 
to carry out for any n > 2 the Farey Series of order n. Using Proposition 12 . 61 we 
will give a second algorithm to compute the Farey series up to any given order n, 
using the continued fraction development. 

Proposition 2.5 (Construction of J- n ). Let n > 2, then the elements of T n , 
{pi/li)i<i<<f>(n)> are recursively defined by: 



(2.16) 



I pi+i = -pi-i + riPi 
I q i+1 = -qi-i + nqi , 



where n = [(n + qi-\)/qi\, starting with = (0,1), (^2,92) = (!,"-) an d 

(P3,?3) = (1,71- 1). 

Proof. Let p/q e T n . Because p and q are relatively prime we can always solve in Z 2 
the linear Diophantine equation qP — pQ = 1: let (Pq,Qq) be a particular solution 
and let r be the integer such that: n — q < Qq + rq < n, namely r = \_{n ~ Qq) / q\ . 

Let us define P r = P + rp and Q r — Q + rq, then the following claims are 
trivial: (P r , Q r ) is again a solution of the linear diophantine equation, {P r , Q r ) = 1 
and < Q r < n. So P r /Q r € T n . Clearly P r /Q r > p/q and we claim that it is the 
immediate successor of p/q in T n . 

To obtain a constructive algorithm we must solve the linear diophantine equation, 
this is achieved by considering the element which precedes p/q in T n : let us denote 
it by p'/q'- A particular solution is then given by Pq = —p', Qo = —q': from 



This is different from the Farey Tree which is still a set of rational numbers in [0, 1] which 
can be constructed by induction starting with: ^0 = {0, 1} and then defining the i— th element of 
f n , n > 1, by: 

The Farey Tree of order n is clearly larger than the corresponding Farey Series and cardj^n = 
2 n + 1. 
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the previous result the element following p/q is then given by: P r = —p' + rp, 
Qr = ~4 + rq, where r = [(n + q')/q\. 

To finish the algorithm we need two starting elements of T n apart of 0/1, but 
it is easy to realize that the first three elements of J- n are 0/1, l/n and l/(n — 1), 
whenever n > 2. □ 

We are now able to give a second algorithm to construct the Farey Series of 
order n. Here is the idea: given an irreducible fraction p/q £ (0,1), we compute 
its continued fraction development and then following two rules: Truncate and 
Subtract one, we obtain two new irreducible fractions in [0, 1] which will result to 
be the predecessor and the successor of p/q in T n with n = q. 

Proposition 2.6 (Construction of JF„, 2 nd version). Let p/q £ (0, 1) and letp'/q' < 
p/q < p" /q" be three successive elements of T q . Assume p/q — [(ciq, Eq), . . . , (an, £«)] 
for some n > 1 and let us define the rational numbers px / qr andps/qs as follows 9 : 

Pt 

(2.17) — = [(a ,£o), ■ • • , (Ofi-i, £s-i)J (Truncate), 

qr 



d 



PS 

(2.18) — = (oo,eo), ■ • ■ , (fln - l,£n)] (Subtract one) . 

qs 



Bn-i = +1. we have 




Pr/qr — p'/q' , and, ps/qs = p"/q" if n is even 
Pr/qr = p" /q" , and, ps/qs = p' '/<?' if n is odd. 
• £n— l = ~ 1 we have the symmetric case, namely 

Pt / qT — p" / q" j and, ps/qs = p'/q' if fi is even 
Pt/ qT = p' / q' , and, ps/qs = p"/q" if n is odd. 



Proof. By fTTQ . (|2~T7|) and l(2~TH|l we have 

\PT = CLn-lPn-2 + £ft-2Pn-3 , | PS = (dft ~ + £ft-lPn-2 

\q T = a,- 1 _iq fi _ 2 + £ft-2<7n-3 1 QS = (an - l)fe-i + £fi-ife-2 



then: 

PT +PS _ an-lPn-2 + £ft-2Pft-3 + (an ~ l)Pn-l + £ft-lPft-2 Pn _ P 

<7t + 9s a ft -i(7 ft _2 + £ft-2<7ri-3 + («s - l)?a-i + Sn-iQn-2 qn q ' 

where we used the definition of p/q with its finite continued fraction of length n. 
Finally: 

P PT _ Pnqn-i -Pn-iqn _ (-l) n+1 £p ■ ■■£„-! 
q qT qnqn-i qnqn~i 

and similarly 

P Ps_ _ Pn(qn ~ qn-i) ~ (Pn ~ Pn-i)qn _ (-l)"£p ■ ■ ■ £„-l 

q qs qnqs qnqs 

from which the proof follows easily. □ 



9 If a n = 2, then En-\ = +1 by remark^] and VsllS = [(oo,eo), ("s-i + 1, +!)]■ 
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2.3. The 1/2— Monoid. In this paragraph we will study the monoid M. oiGL(2, Z), 
introduced in (|2.8(l and used in the construction of the 1/2-Complex Bruno func- 
tion. Our aim is to show its relation with the nearest integer continued fraction: 
for this reason we call it 1/2-Monoid. We will prove that given p/q e [0, 1) we 
can "fill" the matrix = (qlq) in exactly two ways, such that it belongs to M. 
"following the nearest integer continued fraction development" . 

Proposition 2.7. Let p/q G [0, 1), n > 1 and assume p/q = [(a , So), ■ ■ • , (a fi , En)] 
to be the finite continued fraction of p/q. We claim that the matrices gr = (g? g) 
and gs — (qs q)' where the rational pr / qr and ps/qs have been defined in Propo- 
sition \£l\ are given by: 

(2.21) g T = eofl(ai) • ■ ■ g{an-i)£n-ig{an) (Type T) 

(2.22) gs = iog(ai)...g(a n - 1 )in^ 1 g(a n -l)g(l) (Type S) , 

where for i = 0, . . . , n — 1, matrices ii and integer di are defined by: 



(2.23) (oi,£i) = 



(at,l) ifsi = +l 

(a»-l,cr) ifsi = -l 



Before to prove the Proposition we make the following: 

Remark 2.8. It results di > 2 for all i, in fact whenever £j = — 1 one has at > 3 
(see Remark \2.ty) . Because p/q £ [0,1) the first couple (ao,£o) can only be one of 
the following two: (0, +1) if p/q G [0,1/2] or (1,-1) if p/q € (1/2,1). 

Proof. Let k < h and let us introduce matrices £o, . . . , £fc_i and integers So, . . . ,dk 
as in (|2.23J) according to the continued fraction of p/q. Then we claim that g = 
e g(di) . . . g(a k -i)i k -ig{a k ) is equal to {%lz\ P qt ) where p k /q k = [(a ,£ ), ■ ■ ■ , (a k ,e k )}. 
This can be proved by induction (use Remark 12. 81 to prove the basis of induction) 
and then 12.2111 follows by putting fc = n. To prove l|2.22|l it is enough to calculate: 



i Q g(ai) . . .g(a k -i)ik-ig(a k -l)g(l) = 



Pk-l Pk—Pk-l\f0 l\_(Pk-Pk-l Pk 
qk-l q k -q k -i)\l I J \q k -q k -i q k 

□ 



Remark 2.9. Clearly matrices of type T belong to M. (because a« > 2) , whereas 
those of type S belong to the monoid if and only if the continued fraction of the 
rational p/q ends with a couple (an,£n) = (2,1), in fact in this way the matrix g$ 
is given by e g(ai) . . . g(an-i)g(l)g(l) = i a g{di) . . . g(an-i)ag{2), where we used 
that Sfi—i — 1 (because — 2) and ag(m) = g(m — l)g(l) for all m > 2. 

Remark also that, if g is of type T then it cannot ends with ag(2), in fact this 
will implies a continued fraction ending with [. . . , (a ft _i, —1), (2, 1)], but we know 
that this is impossible and so: or £n_i = — 1, &n > 3 and gx = ■ ■ ■ crg(an), either 
e n -i = +1, a n >2 and g T = . . . g(a fi - 1 )g(an) ■ 

With the following proposition we will prove that M. is the union of matrices of 
type T and of type S with (a fi , e R ) = (2, 1). Let us denote by Mt the monoid of 
matrices of type T and Ms those of type S, with (an,£n) — (2, 1). 

Proposition 2.10 (The 1/2-Monoid). M = M T U M s . 
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Proof. Clearly A4t U A4s C M.. Let us prove the other inclusion. Let r > 1, 
mi, . . . , to,- > 2, eo, ... j e r -i € {1, f}, such that g = £o<7(mi) ■ ■ • s r -ig(m r ) € .M. 

Let us consider two cases, first one: £ r -i = a and m r > 3 or e r _i = 1 and 
m r > 2; second case: e r _i = cr and m r = 2. In the former case to g we associate a 
continued fraction by introducing, for i = 1, . . . , r — 1: 

, , f(0,+l) ife< = l , . fK,+l) if e< = 1 

[(1,-1) u £i = o- [(to, + 1,-1) if Si = a 

[(ao,£o), ■ ■ ■ ! (a, r ,£ r )] represents some rational p/q, let us define as before pr/qr 
and then gx = (gj q) = iog(ai) . . .£ r _ip(a r ) (by Proposition 12.71 where we also 
defined a^'s). Observe that Sj = TOj to conclude g = gr ^ Mt- 

The second case can be treated similarly. Now to g we associate the continued 
fraction [(do, £o), ■ ■ ■ , (a r -i, 1), (2, 1)] where, for i = 1, . . . , i — 2: 

f(0,+l) if£i = l , . f(m<,+l) ife t = l 

(a ,£ ) = <, . , (a i ,e l ) = < , a r _i = TO r _i. 

[(1,-1) if £, = o- [(to, + 1,-1) if Si = a 

Let [(ao,£o)i • •• j (or— l) 1); (2, 1)] be some rational p/q, define as before ps/qs = 
[(do, Sq), . . . , (a r _i, 1), (1, 1)], then by Proposition 12. 71 

9s = £os(ai) ■ ■ ■ 9(a r -i)g{l)g(l) = eoff(ai) ■ ■ ■ g(a r - 1 )ag{2) = g 
and it belongs to .Ms- D 

To end this section we introduce a third characterization of the 1/2-Monoid, 
which corrects a small error in §A.4.4 page 836 of [Mi 

rmi et al. 2001| . and which 

will be useful to construct the numerical algorithm for the 1/2-complex Bruno 
function. 

Proposition 2.11. Let g = (" ^) € G. Then g belongs to M if and only if 
d>b>0, c > a > andd>Gc, where Q = (y/E + l)/2. 

The proof can be done by direct computation and we omit it. We end this part 
with the following: 

Remark 2.12 (The Gauss Monoid). In [Marmi et al. 2001] authors considered the 
complex Bruno function constructed using the Monoid M. : 



M 



{•9 = (cd) eG : d>b>a>0 andd>c>a\ 



We recall that according to the Gauss continued fraction algorithm we always have 
ei = +1, we can then prove modified version of Provositions UTo] and \Kl\ to conclude 
that M. is constructed "following" the Gauss continued fraction algorithm: starting 
from p/q £ (0,1), we complete the matrix g* — {qlq) into gs andgx, where ps /qs 
and pt / qr are obtained with the Truncate and Subtract operations acting on the 
Gauss finite continued fraction of p/q. 

2.4. An Algorithm for the 1/2— complex Bruno function. Using the results 
of the previous sections we are now able to give an algorithm to compute the 1/2- 
Complcx Bruno function. Let us rewrite definition (|2.1U|) : 



B(*) = E 



<px/ 2 (z - n) , 
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where ^1/2(2) = — - Li 2 (j-) +i log21og (l — j-) and the action L g has been de- 
fined in (|2.6ll . From the previous sections we know that the sum over M. can be 
replaced by a sum over p/q £ [0, 1) , (p,q) = 1, in such a way that to each p/q we 
associate the matrix g?, and also g$ whenever the continued fraction of p/q ends 
with (a n ,e n ) = (2, +1). 

Using the periodicity and the parity properties of B we can restrict to 5Rz £ 
[0, 1/2]. Let us consider the contribution of some p/q £ [0, 1) to B- Because of the 
form of (pi/2 and of the action L g we remark that the larger is the denominator of 
the fraction, the smaller is its contribution to the sum; moreover different rational 
numbers with the same denominator give comparable contributions, so we decide 
to order the rationals w.r.t. increasing denominators: in other words according 
to the Farey Series. A similar statement holds w.r.t. the sum over Z: large n's 
give small contribution to the sum. We then introduce two cut-off to effectively 
compute H2.10fl : N max denoting the largest order of the Farey Series considered and 
ki the largest (in modulus) n £ Z which contributes to the sum over integers 10 . 




/ // 



/ 



0.1 0.2 03 0.4 0.5 





0.4 0.5 



0.1 0.2 0.3 0.4 0.5 



Figure 1. Plot of B{z) vs dtz at fixed. The top line contains 
SJiB whereas on the bottom line we plot SB- The column a) is 
for = 1CP 3 whereas column b) is for 9z = 10~ 4 . Each plot 
has 10000 points 5ftz uniformly distributed in [0,1/2]. k\ = 80, 
fc 2 = 20, N max = 151. 

Then the 1/2-complex Bruno function can be numerically approximated by: 

/ 

(2.24) B(«)~ J2 J2 L^ : p^{\ + L a )ip l/2 {z-n), 

\n\<ki p/qEFNwax 

where p»/g* £ {pt/qt^Ps/qs} and the sum is restricted to fractions such that 
Q > QlSi (being always q > Gqr)- This approximation can be made as precise as 



^For technical reasons, we prefer to introduce a third cut— off, &2- We refer to Appendix 1X1 to 
explain the role of this cut-off. 
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we want, by choosing N max and k\ large enough, in fact 1(2.10(1 can be obtained as 
double limit N max — » +00 and k\ — > +00. In Appendix lAl we will give numerical 
results showing the convergence of 1(2.24(1 varying the cut-off values, the convergence 
of 3B(z) to B{lRz) when 9z — > and SRz £ B, and the 7r/<7-jumps of 9?B(z) when 
z — > p/g, as proved in [Marmi et al. 2001| . In Figure □ we show some plots of B(z) 
for fixed (small) values of and S [0, 1/2], whereas in Figure El we show two 
polar plots of e zB ^ z ^. 




" fti 005 O; 0.J5 0.2 005 0J 0J5 0.2 



Figure 2. Polar plot of e iB ^^ for fixed values of 3z. a) 3z = 10 2 
and b) 3z = 10~ 3 . 12000 points uniformly distributed in [0,1], 
fcl = 80, fc 2 = 20, Nmax = 151. 

3. The Yoccoz Function 

The aim of this section is to briefly introduce the algorithm used to compute the 
Yoccoz Function, U(X), introduced in § 11.11 Let A e D*, let P\{z) — Xz (1 — z) be 
the quadratic polynomial and let us introduce the polynomials: U n (X) = X~ n P^ n (1/2). 
Then we recall that the Yoccoz function is the uniform limit, over compact subsets 
of D, of U n (X). 

From 11.3(1 and its original definition, H\ (/7(A)) = 1/2, we get: 
(3.1) U{\) = \- n H^ (X n U n (X)) , 

for all integer n. Hence to compute U(X) we need to know how close is H^ 1 
to the identity, near zero and this can be done using some standard distortion 
estimates Bu ff et al. 2001) . So for any fixed A G ©*, we can find n = n(X) s.t. 
P° n (l/2) is contained in some fixed disk on which we can apply the distortion 
estimate and then from (|3.1|l compute an approximation to 17(A) with a prescribed 
precision ejj- 

Remark 3.1 (Parity of Yoccoz's function). Let us observe the following facts. 
Assume X = e 27Tl ( x + lt ) f with t > fixed, and x varying in (0,1/2) and let us 
introduce u(x) — JJ (e 2j ' llyX+lt ^ , to stress the dependence on x only. Then we claim 
that $iu(—x) = $lu(x) and Syu(— x) — —Qu(x). The proof can be done as follows. 
First remark that X, as a function of x, is mapped into X, when x i— ► —x; then 
is enough to observe that polynomials U n (A) verify, for n > 2, U n (A) = U n (A), 
namely: 

W n (e 2vi ( x+it ^ =mj n {e 2 ^- x+it ^ and $U n (e 2wi ( x+it ^ = -%U n (e 2wi< -- x+it A . 
A similar statement holds for \ogU(X). 
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Using the 1,-periodicity , we consider the Fourier series o/C/(e 27T4Z ) and using an 
argument similar to the one of B,emark \2.1\ we conclude that all the Fourier coeffi- 
cients are real and zero for negative Fourier modes. Clearly Taylor's coefficients of 
U(X) coincide with Fourier coefficients ofU(e 2mz ). 

Figure |21 shows some polar plots of t/(e 27 " z ), for different values of 3z > 0, 
whereas in Figure^real and imaginary parts of — log [/(e 27 " z ) are given. Compare 
with Figures |2 and ^ 




Figure 3. Polar plot of U(e 2viz ) for fixed values of $sz. a) 3z = 
10~ 2 and b) 3z = 1CT 3 . 12000 points uniformly distributed in 

[o,i], ec/ -io- 3 . 





Figure 4. Plot of - \ogU(e 2mz ) vs 3?z at fixed 3z. On the top 
we plot the imaginary part whereas on the bottom the real part. 
Column a) is for 3z = 10 -3 , whereas on the column b) we show 
5z = 10~ 4 Each plot has 10000 points uniformly distributed in 
[0,1/2], e v = 10~ 3 . 



Let us conclude this section with the following remark. 
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Remark 3.2. In Figure\^we show some polar plots of the "Yoccoz function" used 
in |Buff et al. 20(1 1] (Figure 2, page 4^4) : they don't look like our previous pictures. 
Here is the reason. They take the following quadratic polynomial Q\ (z) = Xz + z 2 , 
which can be conjugate to our choice, P\(z) = Xz (1 — z), using A(z) = — Xz: 

AoP x = Q x oA. 

Let us denote by V(X) the Yoccoz function for the polynomial Q\, then we claim 
that: 

-XU(X) = V(X) , 

which explain completely the relation between Figure and Figure Because 
— 31ogJ7(A) exhibits the same jumps at rationals as the real part of the complex 
Bruno function does, we choose the quadratic polynomial in the form P\. 





Figure 5. Polar plot of V(e 2mz ) for fixed values of 3z. a) 3z = 
10~ 2 and b) 3z = 1CT 3 . 12000 points uniformly distributed in 
[0,1], e v = HT 3 . 



4. The Littlewood-Paley Theory 

The aim of this section is to introduce the basic ideas and results of the Littlewood- 
Paley Theory, for a more complete discussion we refer to jStein 1970llFrazier et al. 1991j 
and also to |De la Llave-Petrov 2002j where authors apply this Theory to study the 
regularity properties of the conjugating function for critical circle maps. In § 15.11 
we will present the numerical implementation of this Theory to study the Holder 
regularity of the function TL and the obtained estimate for the Holder exponent. 

The decay rate of the coefficients of a Trigonometric series, ^ z Cke 27rlkx , does 
not determine whether this series is the Fourier series of some L p function if p ^ 2. 



,2-Kikx 



then 



More precisely given f £ L p , 1 < p < 2 and its Fourier series ^ z 

for "almost every choice of signs ±1", the series X; z (±l)/ fc e 27 " fc:E is not the Fourier 
series of a L p function. This problem has been overcome by Littlewood and Paley 
by "grouping together" trigonometric coefficients in dyadic blocks. Let A > 1, 
(£ /) (x) = f and, for M > 1, let (C M f) (x) = E^-i<|n|<A« /«e 27rmx , be the 
dyadic partial sum of /. Introducing the Littlewood-Paley d-function: d(f)(x) — 

-, 1/2 

J2m>o \£-Mf(x)\ 2 , one can prove |Littlewood-Paley 1931[ IPrazier et aT~ 1991 
that for all 1 < p < +oo there exist positive constants A p and B p such that: 

„<IK/)||p<Bp 
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The Littlewood-Paley Theory is indeed more general, allowing to character- 
ize other functional spaces by property of Fourier coefficients, for instance it ap- 
plies [Frazier et al. 1991] to Sobolev spaces, Hardy spaces, Holder spaces and Besov 
spaces. In the case of Holder regularity one can easily realize that Fourier coeffi- 
cients of an 77-Holder continuous function decay according to fi = 0(\l\~ v ), the 
converse is not true but again the Littlewood-Paley Theory can characterize the 
Holder regularity by the decay rate of the dyadic blocks. 

An important tool in the Theory of Fourier series is the Poisson Kernel: P s (x) = 
E fcez s^e 2 ™ kx , s G [0, 1) and x e T. Let (/ * g)(x) = ft " # be the 

convolution product for 1-periodic functions. Then one can prove the following 
result f [Stein 1970j Lemma 5 or [Krantz 1983) Theorem 15.6) 

Theorem 4.1 (Continuous Littlewood-Paley). Let < r\ < r. r 6 N and f be a 

continuous 1-periodic function. Then f is i] -Holder continuous if and only if there 
exists C > such that for all t > 0: 



where V f {x,t) = (P exp (- 2 7rt) * /)(«)• 



We remark that if the Theorem holds for some r £ N, then the same is true for 
any r\ € N, f*i > r. We call this Theorem Continuous Littlewood-Paley to distin- 
guish it from the following result, which is more close to the original idea of dyadic 
decomposition and we will call it Discrete Littlewood-Paley (see [Krantz 1983j The- 
orem 5.9) 

Theorem 4.2 (Discrete Littlewood-Paley). Let r\ > and let f e C° (j). Then f 
is "q-Holder continuous function if and only if for all A > 1 there exists a positive 
constant C such that for all M 6 N we have: 

\\C M f Woo < C A-^ M . 

One usually take A = 2, and so the name dyadic decomposition, but the result 
is independent of the value of A. In the numerical implementation of this method 
we will use a value A close to 1.25 for computational reasons. 

5. Presentation of numerical results 

This section collects our numerical results about the Marmi-Moussa-Yoccoz 
conjecture that we recall here: 

The analytic function, defined on H + : z 7i (z) = log J7(e 27r42 ) — 
iB(z), extends to a 1/2-Hdlder continuous function on the closure 
ofW+. 



Let us begin with some consideration concerning TL. Remark l2.1l and Remark l3.ll 
imply that TL{z) has even real part for $tz S [0, 1/2] and an odd imaginary part on the 
same domain. Moreover Tl(z) is 1-periodic and its Fourier series has only real co- 
efficients, which correspond to non-negative Fourier modes: TL{z) — J2i>o hie 2vllz . 

In Figure we plot real and imaginary parts of Tl(z) for some fixed small 3z; 
remark that H still has a "structure" but jumps of B and — log U seem to "com- 
pensate" to give a continuous function. The same fact holds for the "bubbles" of 
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2-i ' 1 ' 1 ' 1 ' 1 ■ r. h 1 1 1 1 1 1 1 1 1 1-2 




°0 0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 L 

Sb . . 




Figure 6. Plot of -iB(z) + \ogU(e 2niz ) vs ftz at fixed 3z. On 
the top we show the real part whereas on the bottom the imaginary 
one. The first column a) is for = 10~ 3 and the second one b) 
for Qz = 10~ 4 . Each plot has 10000 points uniformly distributed 
in [0, 1/2]. fci = 80, fc 2 = 20, N max = 151 and e v = 10~ 3 . 

U{e 2vlz ) and e 1 ^ 2 - 1 , Figure0show some polar plots of e n ^ for fixed small 3z > 0, 
there are still some "bubbles" but they are far from (0,0). 

The Holder continuity will be proved in the next paragraph, by giving an estimate 
of the Holder exponent applying the Littlewood-Paley Theory in the Discrete and 
Continuous versions. 




FIGURE 7. Polar plot of e~ iB(z) [/(e 27riz ) for fixed values of 3z. 
a) 3z = 10~ 2 and b) 3z = 10~ 3 . 12000 points uniformly dis- 
tributed in [0, 1], fa = 80, k 2 = 20, N m = 101, e v = 10~ 3 . 



5.1. Numerical Littlewood Paley method. To numerically implement the Littlewood- 
Paley Theory we compute from the numerical values of 7i, a finite but large, number 
of Fourier coefficients. Then to apply Theorem l4.1l we remark that the convolution 
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with the Poisson Kernel and the r— th derivative have the form: 



' l>0 

where we used the previous remark on the Fourier coefficients of TL. 



We numerically compute 



{ P exp(-2nt) * ft 

t and some r > 1, then applying a linear regression over the data 
d_ 

at 



for several small values of 



(5.1) 



log 



Pex P (~2TTt) * U) = C' r - (3 C Lp(r)\ogt , 



we obtain a numerical value for f\Q LP = r — (3cLp{r). 

From a numerical point of view the continuous version of the Littlewood-Paley 
method is better than the discrete one, in fact the former has two parameters 

(r) 

to vary t and r. We can vary r to control if the computed value of r\(j LP stays 
constant or not. Moreover we can compute the l.h.s. of (|5.1(l for many values of t 
and for each one all the known Fourier coefficients are used, whereas in the dyadic 
decomposition to "small" M only "few" Fourier coefficients give their contribution 
and only for "large" M a large number of Fourier coefficients enter. 

In Figure |H1 we report data from 1)5. Q and the corresponding linear regression 

1 1 (v—2) 

values . The estimated values of i] obtained for different r are: ti clp = 0.497 ± 
0.003, = 0.498 ± 0.004 and rf^p = 0.498 ± 0.003 (errors are standard 

deviation errors of linear regression). They agree in the numerical precision and 
this gives a good indication of the validity of the results. There is no reason to prefer 
one value to the other and so we estimate r\cLP = 0.498 ± 0.004: the mean value 
of the interval obtained by the union of the three intervals obtained for r = 2, 3, 4. 




In the Figure we decided to show only few points to have an "intelligible picture", but the 
linear regression are made using hundred of points. 
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We also report the numerical results obtained using the discrete Littlewood- 
Paley Theorem. We fix some A > 1 12 and from the computed Fourier coefficients 
of TL we construct the dyadic partial sums for some large M £ N- Then we use a 
linear regression on the data: 

(5.2) log^ HAm/Hoo = Cdlp — VdlpM , 

to obtain the estimate value of the Holder coefficient: r]oLP — 0.50 ± 0.03 and 
Cdlp = —4.90 ± 0.66. In Figure El we report data from 1|5.2|1 and the linear 
regression applied on "large M" . 



-w 

^ S -15 
-20 



5 10 15 20 25 30 35 

M 

FIGURE 9. The function M h-> \og A ||£mH||oo and the linear fit 
log^ ||£mH||oo = C DLP - rj DL pM. 



5.2. Conclusion. We conclude this paper by summarizing the obtained results. 
We introduced the 1/2-complex Bruno function and the Yoccoz function, both 
with an algorithm to evaluate them numerically. We studied the function Ti(z) = 
—i B(z) + log U (e 27Tlz ) defined on the upper Poincare plane, and we conclude that it 
can be extended to its closure, with a trace 77-Holder continuous. Numerical results 
based on the Littlewood-Paley Theory give us the estimated value for the Holder 
exponent: t]clp — 0.498 with an error of ±0.004. We can then conclude, with a 
good numerical evidence, that the Marmi-Moussa- Yoccoz Conjecture should hold 
with the maximal exponent 1/2. 

Appendix A. Numerical considerations 

The aim of this Appendix is to consider in details some technical parts of our nu- 
merical calculations. We will consider the role of the cut-off and their relations with 
the accuracy of the computations. We will also compare the numerical properties 
of B with the analytical ones proved in |Marmi et al. 2001 . 




The exact value of A is fixed in such a way we can take M sufficiently large to have a good 
asymptotic, even if we have a finite number of Fourier coefficients. 
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A.l. Accuracy of the algorithm for H(z). Let us recall the formula defining 
the 1/2-complex Bruno function: 



B(z) = J2 



^ L g (l + L a ) 



ip 1/2 (z - n) , 



_g£M 

as already observed we need to introduce three cut-off to compute it: N maxi ki and 
k 2 . The first one determine the largest Farey Series involved, namely only fractions 
p/q s.t. p/q e [0, 1), (p, q) = 1 and q < N max will be considered to compute B- 
The other two cut-off: ki > k 2 > 0, are introduced to truncate the sum over %, 
Because the larger is q, the smaller is its contribution to B, to gain CPU times we 
decide to truncate the sum over Z at \n\ < k\ if q is "small", and to \n\ < k 2 if q is 
"large". Results showed in Section are obtained with N max — 151, k\ = 80 and 
k 2 = 20. 

In the rest of this paragraph we will study the dependence of the computed 
Bruno function on the cut-off. Let us fix two cut-off but one, call it generically 
M, and let us stress the dependence of B on it, by setting Bm(^). We are then 
interested in studying the relative error: e re t(z,M) — |Bm (z) — B (z)|/|B (z)|, 
where z is fixed and B (z) is numerically computed with some fixed large cut-off: 
N m ax = 101, ki = 80 and k 2 — 20. Or we can consider e re i (M) the mean value of 
e re i (z, M) for Jfz £ [0, 1/2] and some fixed value of 9z > 0. In Table |A~T1 we report 
values of log 10 e re i (M), for M 6 {N max , k\, k 2 }. Whereas in Figure ITUl we show, 
log 10 e re i (z, M), for M = N max and M = ki and z E {^2 - 1, 2 - £}. Clearly the 
larger are the cut-off the more accurate are the results, but recall that large cut- 
off implies large CPU times; in particular the CPU times increases almost linearly 
w.r.t. k\ and k 2 , but quadratically w.r.t. N„ 



v max • 




Figure 10. Plot of log 10 e re i (z, M) for some "good" z. On the left 
we show log 10 e rel (z, N max ) whereas on the right log 10 e rel (z, fci). 
Circles are for z = 2 — Q + H0~ 7 and squares are for z = \[2 — 1 + 
il0~ 7 . 



To have a full test of our algorithm we try to evaluate the limit, for 3z — > 0, of the 
computed B and compare it with the results proved in |Marmi et al. 2001] : section 
5.2.9 page 816 and Theorem 5.19 page 827. In particular we will be interested in 
studying the rate of convergence of SB (x + it) to Bi/ 2 (x), for t — > when x is 
some "good" number fFigure ITT)) . and the "jump value" of 5i B (p/q + it), when t 
is "small" (Table A.2). 



20 



TIMOTEO CARLETTI 



N 

1 v qmax 


l°gl0 e rel {Nrnax) 




log 10 e re i (fci) 


fc 2 


log 10 e re i (k 2 ) 


83 


-5.90 


60 


-9.57 


15 


-10.80 


61 


-4.73 


40 


-9.10 


10 


-10.32 


41 


-3.81 


20 


-8.63 







Table 1. We report log 10 e rei (N max ), log (fci) and log 10 e re i (k 2 ). 



-.-2 

t-3- . 

d ' 

" 5 -S -7 -6 -5 -4 -3 -2 

FIGURE 11. Convergence of 3B(s + it) to Bi/ 2 (x) when f — * 
for some "good irrational' x. Circles are for x = 2 — C7 whereas 
squares for x = \[2 — 1. We plot log 10 |^B (x + it) — B 1 / 2 (x)| 
versus log 10 t. N qmax = 101, fci = 80 and k 2 — 20. 





A5iB (p/q + it) — ir/q 


0/1 


1.1 10~ 3 


1/2 


7.4 10" 4 


1/3 


3.6 10" 3 


1/4 


3.610" 3 


1/5 


3.8 10" 3 


2/5 


2.5 10" 3 



Table 2 . The jumps of 5ft B (x + it) for rational x and small <. The 
jump at x = p/q is the numerical difference |5ftB {p/q + S + it) — 
5ft B {p/q — 5 + it) |, for 5 small. We report the difference of the 
jump w.r.t the expected value for x G {0/1, 1/2, 1/3, 1/4, 1/5, 2/5} 
and t = 1Q- 7 . 
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